The strong-coupling master equation of quantum open systems 
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In this paper we demonstrate how to generate the strong-coupling master equations for open 
quantum systems of continuous variables. These are the dissipative master equations of quantum 
Brownian particles for which the environmental noise is stronger than other system forces. Our 
strong-coupling master equations are very different from other so-called "strong-coupling" master 
equations (e.g. the quantum Smoluchowski equation) which are perturbing off a limit in which the 
system energy is taken to be perturbative and thus the dynamics is principally Markovian. Such 
approximations also require the system mass to be asymptotically large (even as compared to the 
ratio of noise and induced system frequencies) and thus they do not fully categorize the regime of 
what one might consider to be strong coupling. Our master equations are highly non-Markovian and 
radically different for different system potentials, admitting no apparent generic form. This result is 
quite exciting as it brings forth a new regime for theoretical exploration: the regime of strong noise 
and dissipation yet non-Markovian, such as strong coupling to a low-temperature environment with 
large 1/f fluctuations. 
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I. INTRODUCTION 

An open system is obtained by coarse graining over 
the environmental degrees of freedom, and its dynamics 
is described by a master equation governing the reduced 



density matrix or marginal Wigner function. The coarse- 
grained environment can act as a source of noise, dissipa- 
tion, and decoherence; thus its influence provides a mi- 
croscopic model of dissipative quantum mechanics more 
general than Markovian (white-noise) models which can 
be constructed more phenomenologically. Exact master 
equations for the stochastic dynamics of open quantum 
systems are, in general, beyond the reach of the simple 
theorist. The canonical regime of interest for nonlinear 
quantum open systems has been one of weak coupling be- 
tween the system and environment. The resultant per- 
turbative master equations for the dynamics of the re- 
duced system are well known and can be derived in a va- 
riety of ways [l-'-i]- At zeroth-order the dynamics are that 
of the uncoupled system, and further orders consistently 
introduce noise and backaction from the environment. 

At the other end of the spectrum, there is also great in- 
terest in the strong-damping regime. Naive expectations 
of behavior in this regime appeal to the limit in which 
the system Hamiltonian is taken to be vanishing and the 
system is primarily driven by the system-environment in- 
teraction. For a single (collective) coupling to the envi- 
ronment, or univariate noise, this produces a rapid de- 
coherence in the basis of the system coupling operator. 
Ergo, for a strong position coupling to the envirormient, 
it is thought that the system should decohere into a po- 
sition distribution. However, this limit has questionable 
applicability to small mechanical systems. Large noise 
fluctuations can generate large kinetic energies, which 
are then incompatible with the limit of vanishing system 
energy. 

In this work we derive from first principles the mas- 
ter equation of a quantum system primarily driven by 
noise and dissipation. What we take to be perturbative 
is not the entire system Hamiltonian as with the "strong- 
coupling" quantum Smoluchowski equation[-i, 5] but only 
the system potential, or nonlinearities thereof. The re- 
sult is a class of master equations which are highly non- 
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Markovian and seem to admit no generic form. Instead 
we have a procedure for generating said master equa- 
tions through an iterative appUcation of noise averaging 
and integration. 



II. OPEN SYSTEM REVIEW 
A. Classical-like Brownian trajectories 

Let us consider the canonical non-Unear Brownian mo- 
tion problem: a Brownian particle (the system) with 
phase-space coordinates x, p and potential V{x) bilin- 
early coupled to an environment consisting of many os- 
cillators with coordinates Xk and pk ■ Linear coupling to 
a bath of oscillators provides a general model of Gaus- 
sian noise [6] . Although we take the dynamical equations 
to be of the classical variety, we will choose a quantum 
initial state for the environment. The reason for this will 
become clear when we derive our perturbative formalism. 

We can write the full system + environment Hamilto- 
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with Ck the coupling strength and where the system- 
environment coupling has been inserted in the bath po- 
tential as a method of renormalization: in this manner 
the phenomenological linear spring constant of the sys- 
tem is the same for any system-environment coupling Cfe. 
The system and environment are taken to be initially un- 
correlated and uncoupled, with the environment initially 
being in its free (quantum) thermal state. The interac- 
tion is switched on instantaneously via Ck at the initial 
time. The equations of motion for the system are then 



x{t) = -Pit) 



(11.2) 



p{t) = -V'[x{t)] - mSculxit) + ^ Cfc Xk{t) , (IL3) 



and the equations of motion for the environment are 



ik{t)^ Pk{t), (n.4) 

Pk{t) = -rukojl Xk(t) + Ck x(t) . (IL5) 

For each equation the momentum can be substituted by 
position to obtain 

mit{t) + V'[x(t)]+m5uj1^x{t) = ^CfcXfc(i), (IL6) 



ruk Xk{t) + rukujl Xk{t) = Ck x(t) . 



(11.7) 



Additionally, the environment positions can be solved in 
terms of the system position 

Xk{t) = nikMt) Xk{0) + fk{t)Pk{0) + Ckfkit) * x{t) , 

(11.8) 
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where * denotes the Laplace convolution A{t) * B{t) — 
dr A{t—T) B{t) and fk {t) is the free Green's function 
of an environment oscillator. This expression can be sub- 
stituted back into the equation of motion for the system 
position, providing the quantum Langevin equation 

miit) + V'[x[t)] + 2m7(t) * x{t) -f 2m'^{t) x{Q) = ^(i) , 

(11.10) 

where the positive-definite, stationary damping kernel 
7(i, r) — j{t — T) and Gaussian noise process £,{t) are 
respectively given by 



2m ^ — ^ mkCdf 

k ^ 

at) = Y.Ck [mkMt)xkiO)+fkit)pkiO)] . (IL12) 



The stationary noise kernel is defined 

Ht,T) = mar))^ = H~r) 



(IL13) 



and is related to damping by the (quantum) fluctuation- 
dissipation relation 



^{e) = m 7(e) e coth^ 



2TJ 



(11.14) 



where 7(e) = dte~^'^* ■^{t). Essentially, the damping 
kernel (equivalently the spectral-density function) and 
temperature completely characterize Gaussian, thermal 
noise. 



1. Linear systems 

For a linear system the Langevin equation takes the 
classical form 

m x{t) + muj^ x{t) + 2mj{t) * x{t) + 2mj{t) x{0) ^t) , 

(11.15) 

and the solutions are given by 

x{t) = mgit) xiO) + git)piO) + git) * ^(0 , (11.16) 

where the Green's function g{t) is most easily determined 
in the Laplace domain as 



s2 + Cj2 + 2s7(s) ' 



(11.17) 
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where g{s) = dte~^* g{t). Defining the phase-space 
coordinates z and homogeneous system propagator $(t) 
as 

zit) = [xit),p{t)f , (11.18) 

T% 1 , (11.19) 

m'^gyt) mg(t) ^ ' 



with free environment propagator 

mkfk{t) fk{t) 
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(11.20) 



we can express the system phase-space solutions 



z(<) = *(t) z(0) + *(t) * ^^kit) Zfc(O), (11.21) 

k 

(11.22) 
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and the homogeneous system propagator $(t) can be 
identified as the system-system portion of the full system- 
environment propagator. Similarly the environment so- 
lutions can be expressed 



Zfe(i) =cl),(0*c,*(t)z(0) 



(11.23) 
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[5kA^ + *fe(i) * Cfe *(i) c,] * </,^.(t) z,(0) . 



The stochastic description extends to our phase-space 
representation with stochastic process H(t) = [0,^(<)]"'" 
in phase-space coordinates, such that 

z(t) = *(t)z(0) + *(i)*H(t), (11.24) 

where the Gaussian noise correlation is now given by 
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where v{t,T) is the noise correlation of ^(t) given by 
Eq. (11.14). 



where the homogeneous propagator #(i) is given by 
Eq. (11.19) and the thermal covariance is 

axit) = d.Ti dT2 *(i - Ti) v{ti,T2) [t - ra) . 
Jo Jo 

(11.30) 

This is a special case of the two-time thermal covariance 

(TT{tl,t2)= / dri / (ir2*(ti-Ti)l/(Ti,T2)*'^(i2-T2), 

Jo Jo 

(11.31) 

which is the thermal contribution to the two-time corre- 
lation function 

(z(<i) zT(t2))^ ^ $(fi) o-o *^(t2) + MtiM) , (11.32) 

given the initial correlation ctq = (z(0) z'^(O)). The first 
contribution is thus homogeneous and dissipative. 



III. A NEW PERTURBATIVE FORMALISM 

A. Classical Dynamics 

Let z = {x,p) and z^. — {xk,Pk) denote the individual 
system and environment phase-space coordinates respec- 
tively. Let Z = (z, zi, Z2, • • • ) denote the collective sys- 
tem -I- environment phase-space coordinates. The classi- 
cal dynamics of phase-space distributions is such that 



W[Z{t),t\ = iy[Z(0),0] 



(IILl) 



where Z(i) is a classical trajectory with initial conditions 
Z(0). Therefore the classical propagator Gc(ti,t2) : 
W{Z,t2) — > W{Z,ti) must act upon arbitrary distribu- 
tions d such that 



B. Linear quantum Brownian motion 

For a linear system, the master equation [ ] of the re- 
duced density matrix p can also be represented in the 
phase-space representation [ ] as a Fokker-Plank equa- 
tion of the reduced Wigner function W{z,,t), given by 

jW{z,t)^Cit)W{z,t), (IL26) 
£(<) = V/ n{t) z + D(0 V, . (IL27) 

The Fokker-Planck equation takes a classical form except 
that the diffusion coefficients describe quantum fluctua- 
tions. The time-local homogeneous and diffusion coeffi- 
cients are given by 

'H{t) = -^{t)^-\t), (IL28) 

D(i) ^]^[n{t)cTT{t)+a-T{t)'H^{t)+&T{t)] , 

(11.29) 



Gc{tiM)d[Z{ti)]^d[Z{t2)]. [111.2) 

The classical trajectories are determined by classical 
equations of motion of the form 

|z(t)=F[Z(t),t], (IIL3) 

and correspondingly, the classical Fokker-Plank equation 
is given by 

^Mz, t) = -F(Z, i)T Vz W{Z, t) . (in.4) 

In contrast, the quantum master equation is not generally 
linear in derivatives, thus not admitting trajectories or 
"characteristics". For polynomial potentials beyond the 
quadratic order, the quantum equation contains higher- 
order derivatives accompanied by corresponding factors 
of h. In general the time-translation generator is nonlocal 
in the phase-space coordinates coordinates ['•]. 
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B. Quantum dynamics in the characteristics 
picture 

We now examine the quantum dynamics of the com- 
bined system + environment generated by a Hamilto- 
nian, Eq. (II. 1), which yields unitary evolution. In the 
quantum phase-space representation, the Wigner func- 
tion obeys 



d_ 
di 



WlZ,t]^Cz[Vz,Z,t]W[Z,t], (111.5) 



where Cz denotes the unitary generator of the system -|- 
environment (compare with Eq. (11.26) of the open sys- 
tem). Let us inspect the dynamics of quantum states 
along some curves Z(i) which are not assumed to be 
proper characteristics. 



dt 



W[Z(t),<] = 



1 T 



dt 



Vz(t)VK[Z(0,t] (III.6) 



£z[Vz(t),Z(i),t] W[Z(i),i] 



If these curves have some first-order equations of motion 
governed by the "forces" F{Z,t) then we can express the 
dynamics along these curves as 



dt 



Wim.t] = 5C[Vz(t),Z{t), t] W[Z[t),t] , (111.7) 



5£(t) =£z(t) + F(Z,t)TVz 



(Ill.i 



The key point is that we will utilize curves Z{t) such 
that 6C{t) is ordinarily a system operator and then cal- 
culate the corresponding perturbative open-system mas- 
ter equation for the marginal distribution W{z,t), which 
is equivalent to the reduced density matrix. 

Any Gaussian influence acting on the system can be 
modeled via linear coupling to a linear environment, 
and the dynamical contributions of all such terms are 
first-order (dynamically classical) in the phase-space rep- 
resentation. Therefore we can always transform to a 
classically-evolving coordinate system, along which the 
dynamics of the environment are effectively integrated 
out. Moreover, the system Hamiltonian used for the clas- 
sical curves is irrelevant for the specific purpose of inte- 
grating out the environment dynamics. We do not have 
to use the true system Hamiltonian for our characteris- 
tics, but for the purposes of perturbation theory we do 
want SC{t) to be small. Therefore we generally want the 
classical dynamics to reproduce as much of the quantum 
dynamics as possible, allowing for much cancellation in 
Eq. (III. 8). Essentially, the noise average of Eq. (III. 7) is 
the nonlinear generalization of the stochastic description 
of QBM found in Ref. [10]. Whereas for linear systems 
the stochastic trajectories z{t) immediately provide the 
solution, here they only serve as an evolving coordinate 
system which encapsulates all effects of the environment. 



1. Perturbation along characteristics 

The integral equation of motion corresponding to 
Eq. (III. 7) is 



W[Z{t),t] =WIZ(0),0] 



(III.9) 



/ dT6C[Vzir),Z{T),T] W[Z{t),t]. 

Jo 



which is amenable to perturbation via a Neumann series 
yielding the perturbative solutions 



Wo[Z{t),t]^W[Z{0),0], 



(III.IO) 



Wi[Z(t),i]= / dr,5£[Vz(.),Z(T),r]W[Z(0),0]. 
Jo 

(III.ll) 

and so on, where W[Z,t] = S^o^"!^'*] ^'^^ 
Wn[Z,t] — 0((5£"). Then we can apply the classical 
propagator of Eq. (111.2) to transform back to the (ini- 
tial) domain coordinates. 



VKjZ(0),i] = Gc(t,0)W„[Z(0,i] 



(III.12) 



Next we assume the initial state of the system -I- environ- 
ment to be an uncorrelated product of marginal distribu- 
tions, W{Z, 0) = W{z, 0) Hfc W{zk,0). Upon tracing over 
the environment we obtain the perturbative open-system 
propagator 



Go(t) = (Gc(i,0)) 



E ' 



Gi(i) = fdT (Gc(t, r) SCit) Gc(t, 0))^ 
Jo 

and so on for G„(t) such that 

oo 

G{t) = Gn{t) : W{z, 0) ^ W{z, t) , 

n=0 

and where the operator 

SC{t) = SC[Vz.Z,t] = 5C[V^,z,t], 



(III.13) 
(111.14) 

(III.15) 
(III.16) 



is always taken to be exclusively a system operation. Per- 
turbative solutions of this form are inherently secular in 
time and do not respect Lie group symmetries such as 
unitary or completely-positive semi-group evolution. In- 
stead we will only use this expansion to obtain the per- 
turbative generators. Using ordinary perturbation the- 
ory, our open-system master equation is then 



Co{t) 



|(Gc(M)). 



(Gc(i,0))-\ (IIL17) 



CAt) = SC{t)+jy!^j^ - Ad[Coit)]^ SC{r,t) , 

(IIL18) 

M(T,t) = (Gc(t,r)<5£(r)Gc(T,0))E(Gc(i,0))-^ , 

(III.19) 
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and so on, where C{t) = X]^o ^ri{t) is the open-system 
Liouvilhan or time-evolution generator, such that 



d_ 
di 



Wlz,t] ^C{t)W[z,t]. 



(III.20) 



The lowest-order propagator Go(i), given by 
Eq. (III. 13), is the evolution operator for the clas- 
sical Brownian motion problem, but with quantum 
noise, making it potentially exact for linear systems. 
Co{t), given by Eq. (III. 17), is the corresponding Fokker- 
Plank equation but with quantum noise, e.g. the HPZ 
equation [7]. Essentially, if one can solve the classical 
Brownian motion problem (with quantum fluctuations) 
then one can determine the quantum corrections due to 
the combination of nonlinearities and noise. If one is in 
the classical, linear, Markovian or noiseless regime, then 
all perturbative corrections can vanish in this formalism 
given the appropriate characteristics to perturb from. 

Our hierarchy of approximation schemes is therefore 
not only determined by the order to which we calculate 
the master equation, the n of Cn , but also the amount of 
detail about our system captured in the characteristics 
we use. Perturbing off the damped (but otherwise) free 
characteristics is strong-damping perturbation as JCfji is 
0{V'^). Perturbing off of the damped oscillator charac- 
teristics is a kind of quasi-linear feedback approximation. 
And perturbing off of the full classical characteristics is 
a kind of semi-classical feedback approximation. At least 
formally, one might also imagine an extension beyond the 
classical characteristics into the quantum characteristics 
[11], possibly with an additional (partial) semi-classical 
expansion therein. Each class of characteristics encom- 
passes the previous. 

Note that the system operation dC{t) is fully retained 
in the first-order correction Ci{t); only its effect through 
the environment (back-action or feedback) is approxi- 
mate. To zcroth order, the environment perceives the 
system evolution as being linear or classical and its feed- 
back upon the environment is therefore approximate in 
that manner. At first order we pick up the full nonlin- 
earity of the system in SC{t) along with feedback cor- 
rections from the time integral in Eq. (III. 18). Essen- 
tially SC{t, t) is a two-time and open-system generaliza- 
tion of SC{t), such that ^(t,t) = dC{t) and SCiO,t) = 
Go{t) 5C{t) G'(^^{t), which is integrated over the past in 
order to compensate for what is missed in the lower-order 
approximation. 

This is an analogous formalism to the more common 
weak-damping perturbation wherein the system and en- 
vironment are uncoupled at zeroth order. In weak- 
damping perturbation the higher-order generators in- 
volve integrals over the past which contain the free prop- 
agators acting upon the system- environment interaction. 
Here in strong-damping perturbation the higher-order 
generators involve integrals over the past which con- 
tain the approximate (yet interacting) propagators acting 
upon the system nonlinearities. 



IV. LINEAR BACK-ACTION 

Perturbing off a linear system -I- environment is the 
most straightforward application of this formalism. This 
includes both the quasi-linear feedback and strong- 
coupling regimes. In evaluating the first-order Liouville 
operator in Eq. (III. 18), we aim to calculate the two-time 
and open-system operator 

SC{T,t) ^ {Gc{t,T)SC{T)Gc{r,0))^{Gc{t,0))^' , 

(IV.l) 

for all ordinary system operators 6C{t). This turns out 
to be a complicated procedure and therefore there will be 
no simple strong-coupling master equation which we can 
write down for an arbitrary system potential. Instead one 
must apply the following procedures iteratively and work 
out the master equation for specific system potentials. 



A. Evaluation of two-time open-system operators 

1. Transformation of Derivatives 

As a simple example, let us consider the reduced (lin- 
ear) derivative operator 



5C 



(IV.2) 



SC{r,t) = (Gc(t,T) VjGc(r,0))^(Gc(i,0))E' ■ 

(IV.3) 

First we move all system derivatives to the left by per- 
forming the transformation 



(IV.4) 



Gc(i,T)VjGc(T,i)l Gc(t,0)\ (Gc(<,0))e^ 



System derivatives transform with the transpose of the 
system -I- environment propagator 



V/ ^ Vj $(t - r) V [* * Cfc cljfe] {t - T) 



system 



k 



environment 



(IV.5) 



where we have used the system portion of the full system 
+ environment propagation detailed in Sec. II A 1. Ho- 
mogeneous system derivatives pass through the noise av- 
erage, whereas noise derivatives will be shown to trivially 
vanish in the noise average. Integrating over environment 
derivatives results in only boundary terms which should 
limit to zero. E.g. 



^ dxk^W{xk) = W{xk)\ 

-oo OXk 



+00 



= 0. 



Therefore our reduced derivative is simply 
when left of all other operations. 



(IV.6) 



(IV. 7) 
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2. Transformation of coordinates 



3. Noise averages 



As a simple example, let us consider the reduced (lin- 
ear) coordinate operator 

dC = z., (IV.8) 
^(r,i) = (Gc(i,T)zGc(r,0))E(Gc(i,0))-^ . (IV.9) 

First we move our system coordinates to the right by 
considering the transformation 

6C{T,t)= (IV.IO) 
{Gc{t, 0) [Gc(0, r) z Gc(t, 0)])^ (Gc(i, 0))'^ . 

System coordinates transform with the inverse of the 
propagator (11.24) and so we have 

SC{T,t)= (IV. 11) 

(Gc(i, 0) [*(t) z + *(r) * H(r)])3 (Gc(t, 0))^^ . 

in the stochastic representation. Homogeneous system 
coordinates can pass through the noise average, but then 
exist sandwiched between open-system propagators. To 
simplify the homogeneous part, we construct the generic 
identity 

Az = Az(Gc(t,0))3(Gc(t,0))3^ , (IV.12) 

which evaluates to 

(Gc(t, 0) [A *(t) z + A *(t) * S{t))^] (Gc(i, O))^^ , 

(IV.13) 

and then we choose A = $(r, t) so that the homogeneous 
part of Eq. (IV.13) matches with the homogeneous part 
of Eq. (IV.ll). Therefore we now have 

SC{T,t) ^ ^{T,t)z (IV. 14) 

- (Gc(t,0)*(T,t)*(i) * H(t))3 (Gc(i,0))-^ 
+ (Gc(t, 0) *(t) * H(t) )3 {Gc{t, 0))^' . 

This expression can also be written 

^(T,i) = *(r,t)z (IV.15) 
'Ccit, 0) ^ dr' *f (t, t') H(r')y (Gc(i, 0))^^ , 

where the $f(T, r') is the final- value propagator [o], and 
is given by 

*f (t, t') = e{T-T') *(r-r') - *(r, i) *(i-T') . 

(IV.16) 

This would be the advanced propagator for local dissipa- 
tion 

*adv(T,T') = -0(r'-r)*(r-r'), (IV.17) 
but for nonlocal dissipation no such propagator exists. 



Tracing over environment coordinates will yield mo- 
ments of the noise. Most easily, the environment co- 
ordinates can be identified with the stochastic process 
H(t), a Gaussian noise process with noise correlation 
{S,{t)S^{T))s = i^{t,T) given by Eq. (11.25). There- 
fore we can evaluate our noise averages with the help of 
Novikov's formula (functional integration by parts), e.g. 



(Gc(t,0)H(r'))= = 
- J^dr" /Gcit,0)uiT\T") 



(IV.18) 



Sz{t) 



^z(t) 



_(5H(t")_ 

where from Eq. (11.24), the functional derivative must be 

= *(f-T"), (IV.19) 



Sz{t) 



_SS{t")_ 

and so for calculating (IV.18) we have 

(Gc(t,0)H(r'))3= (IV.20) 

- /dT'V(T',T")*^(i-r")V.(Gc(t,0))3 , 
Jo 

upon substitution and transforming the derivative 
through the propagator. All together the co-rotating co- 
ordinate becomes 



K{T,t) = *(r,i)z 



(IV.21) 



- f dr'l dr"*f(T,T')i^(r',r")*'^(t-r")Vz, 

JQ Jo 

which can be simplified to 

SC{t, t) = *(r, t) z + [*(t, t) (TT{t, t) - ,tt(t, t)] , 

(IV.22) 

where (Tx{T,t) is the two-time thermal covariance in 
Eq. (11.31). As our noise is Gaussian, all operations can 
be reduced to integrals over the homogeneous propagator 
and two-time thermal covariance. 

B. Consistent results 

1. External Forcing of an Oscillator 

Consider the quantum damped oscillator with exact 
master equation given by Eq. (11.27) under the pertur- 
bation of an external force 



(IV.23) 



in the full unitary theory. It was proven in Ref. [b] 
that this is not the force one finds in the open-system 
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master equation. In evaluating the first-order correc- 
tion, Eq. (III. 18), we must invoke our rules for trans- 
forming system derivatives. Here we can specifically use 
Eq. (IV. 7) to obtain 



dC(r,t)= - Vj*(i-T)F(r), 
Ci{t) =6C{t) 



(IV.24) 
(IV.25) 



- j^T I ^ - Ad[£o(t)] I V,^ *(t-r) F(r) . 

The diffusion generator in CQ{t) (11.27) commutes with 
the external forcing, but the homogeneous generator does 
not commute and one easily finds 

£i(t) = 5C{t) - V/ f^dr 1 1 -f -Hit)^ *(t-r) F(r) . 

(IV.26) 

which is exactly the correct effective force; note that this 
correction vanishes for local dissipation where H{t) is 
constant and ^{t) ~ ~H'^{t). The true force is not 
seen in the non-Markovian master equation because, de- 
spite its Markovian representation, the response is still 
inherently nonlocal. 

2. Linear forcing of a free particle 

Finally we consider a peturbation of linear forces which 
take the form 



<5;£(t) = VjK(t)z, 



(IV.27) 



so that the perturbative solutions (11.20) to the Langevin 
equation are simply given by 

* = *o + *i + • • • , (IV.28) 
*i = - / dr *o(i-r) K(r) *o(t) . (IV.29) 



This is sufficient to calculate the perturbative master 
equation using the formulas II B. The simplest coeffi- 
cients to calculate are the homogeneous coefficients 



n{t) = nait) + Tiiit) 

?ii(t) =K(t) 

d 



dr 



dt 



noit) 



(IV.30) 
(IV.31) 

*(t-T)f>K(r)*(T,t), 



which would be trivial in the limit of local dissipation. 
Our new perturbative formalism requires evaluation of 
the reduced operator 



Gc{t, t) V/ K(t) z Gc(r, 0))^ (Gc(i, 0))^' , 

(IV.32) 



which from the results of Sec. IV can be shown to be 

dC{T, t) = Vj *(t - r) K(t) (IV.33) 
X {*(r, t) z + [*(r, t) (TT{t, t) - (Tt{t, t)] VJ . 



It is then a straightforward calculation of Eq. (III. 18) to 
show that the homogeneous terms here exactly reproduce 
those generated by (IV.29). Comparison of the diffusion 
coefficients is considerably more taxing, but they also 
work out to be exactly the same. 



C. New results 

From Sec. IV B we now have for the external driving 
force 

8C{t) = -V^nt), (IV.34) 
^(r,t) = -Vj*(t-T)F(t), (IV.35) 

whereas for the linear force 

J£(t) = V/Ko(t)z, (IV.36) 

we have the two-time operator (IV.33). In any case, the 
first-order master equation perturbative in these forces is 
given by 



£i(<) =^£(t) + ^dr 1^ - Ad[£o(i)]|^ 



(IV.37) 



d_ 
dt 



Ad[Co{t)]^ j^dTSC{T,t), (IV.38) 



in terms of their respective two-time open-system forces. 
These results are compatible with known results for the 
linear system and now we will proceed to derive new re- 
sults. First we consider an arbitrary quantum deforma- 
tion 

5C{t)^^^Fdb{t)x\ (IV.39) 
and then express it in matrix notation as 

6C{t) - (v/p)' F^,{t) (x^zf , (IV.40) 

where x and p are unit vectors in phase space. Applying 
the rules of Sec. IV A, the two-time operator takes the 
form 



8C{r,t)^ [Vj*(i-T)p 



k=0 



(r, t) 



Fdbir) (IV.41) 
Y^Q [i-*(r,OzrAW(r,t), 

x^/" dr' *f(r, r')H(T') 
Jo 

^{Gc{t,0))s\ (IV.42) 



^Gc(i,0) 






. Jo 



where we have expanded the transformed coordinates via 
the binomial theorem. The final step is to evaluate the 
propagated noise moments via successive application of 
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Novikov's formula. The zeroth moment is trivial and the 
first moment we have already calculated 

A[°l(T,t) = l, (IV.43) 
aI^I (t, t) - ±^ [*(r, t) rrrit, t) - ctt(t, t)\ V, , (IV.44) 

while all higher-order operators can be determined recur- 
sively from integration by parts and the product rule. 



and in matrix notation 



aI*^! (r, i) = aW (t, Al'^^il (r, 

-(fc-l)s(T,t)A['=-2l(r,t), 



(IV.45) 



s(t, t) = ^ \ctt{t, t) + *(t, t) arit, t) *^(t, t) 



-2^[^{T,t)(TT{t,T)]±, 



(IV.46) 



and so^^\t, t) is a fc -order differential operator. Thus 
we have explicit rules for all master equations with poly- 
nomial potentials. 



1. Quadratic forcing 

A quadratic correction to the spring force (with some 
cubic correction also in mind) can be written 



.£(t) = fc,(t)(|-.^-l|l 



(IV.47) 



where the second term arises only in the quantum Fokker- 
Planck equation; it carries an additional dimensional fac- 
tor of and generates the quantum deformation of tra- 
jectories. To simplify calculation we rewrite these terms 
in matrix notation as 

6C{t) = k,{t) (v/p) (x^z)^ - ^ (v/p)' . 

(IV.48) 

Applying expansion formula (IV. 41), we find the two- 
time operator to take the form 

dC{r.t) = fci(T) (V7*(t-T)p) (x'r*(T,t)z)' 

+ 2fci(T) (Vj*(t-T)p) (xT*(T,t)z) AW(T,t) 

+ k^{T) (v/*(t-r)p) f AW(T,t)2 - s(r,0' 
ki{r) 
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V,^*(t-r)p 



2. Cubic forcing 



(IV.49) 



The cubic correction to the spring force can be written 
^£(^) = A.(^)(|.^-i|^.), (IV.50) 



6C{t) = Ut) (Vjp) (x^z)^ - M!2 (v/p) (x^z) , 

(IV.51) 

which yields the two-time operator 
5C{r,t) - fc2(T) (v/*(<-r)p) (xT*(r,i)z)' 
+ 2 fc2 (t) {y^^{t-T) p) (xT#(r, t) z) ' All! (r, t) 
+ 2fc2(r) (Vj*(t-r)p) (xT*(r,i)z) (a^ (r, t)^ _ s(r, i) 
+ hir) (v/*(t-T) p) ( aIH (r, 0^-3 s(r, t) aI^ (r, t) 



Vj*(t-r)p) (iT*(r,t)z) 



fc2(T) /_T 



(v/*(<-r)p) AW(T,t). 



(IV.52) 



V. DISCUSSION 



We have derived a fairly general theory of strong- 
coupling perturbation for continuous variable systems. 
Our formalism makes heavy use of the phase-space repre- 
sentation as the zeroth-order problem is a linear one and 
the phase-space representation is perhaps the simplest 
formalism for that case. Essentially we are perturbing 
off of QBM, where the QBM-like model provides us with 
stochastic coordinates which do not solve the problem, 
but integrate the environment dynamics. This results 
in a combined nonlinear and stochastic dynamics of the 
unraveled system. From there we apply the standard 
perturbation theory of master equations, one wherein 
the perturbative time-translation generator is calculated. 
This method retains important Lie-group symmetries 
and grants the possibility of non-secular behavior (thus 
allowing for late-time solutions). As one might expect 
our master equation is highly non-Markovian and highly 
model specific. 

These result are given in a form which is calculable, 
but far more complicated than what one would desire at 
first order. At least two avenues of attack must be con- 
sidered. In the larger view one must question if there is 
a more suitable formalism which includes all necessary 
ingredients: the stochastic map, open and closed-system 
propagation, etc.. The linearity of the zeroth-order dy- 
namics is simplest here in the phase-space representation, 
but perhaps we could trade some of that simplicity for 
a representation which is more apt for the nonlinearity. 
Or perhaps a combined approach could best utilize prop- 
erties of different representations. On the other hand, 
in the more localized view one must address some sim- 
plified calculation of the noise averages which predomi- 
nantly rely upon knowledge of the two-time thermal cor- 
relation and integrals thereof. We have demonstrated 
that a true strong-coupling master equation exists, and 
we have explicitly given its form, but a large amount of 
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work remains to apply our results to representative phys- ical setups. 
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